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Abstract 



^^ Stability analysis and control of linear impulsive systems is addressed in a hybrid framework, through the use of continuous- 
\^ ' time time-varying discontinuous Lyapunov functions. Necessary and sufficient conditions for stability of impulsive systems 
with periodic impulses are first provided in order to set up the main ideas. Extensions to stability of aperiodic systems under 
minimum, maximum and ranged dwell-times are then derived. By exploiting further the particular structure of the stability 
j^ ' conditions, the results are non-conservatively extended to quadratic stability analysis of linear uncertain impulsive systems. 
These stability criteria are, in turn, losslessly extended to stabilization using a particular, yet broad enough, class of state- 
feedback controllers, providing then a convex solution to the open problem of robust dwell-time stabilization of impulsive 
systems using hybrid stability criteria. Relying finally on the representability of sampled-data systems as impulsive systems, 
the problems of robust stability analysis and robust stabilization of periodic and aperiodic uncertain sampled-data systems are 
►^ ' straightforwardly solved using the same ideas. Several examples are discussed in order to show the effectiveness and reduced 
j-y-) , complexity of the proposed approach. 
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1 Introduction 



Impulsive systems [1, [H, \lB, \lB, IHj IH, IH, H, Q are an important class of hybrid systems exhibiting both continuous- 
d • and discrete-time dynamics. The discrete-time part, which is only active at certain time instants t^, fc G N, introduces 
discontinuities in the overall trajectories of th e sy stem. They arise in several fields such as epidemiology [62, [l3| , 
sampled-data and networked control systems [H, HH, [H, 133 1^. [isl , reset systems [H, 0, [13 , forestry ]67| , power 
electronics [4J|, harvesting problems [TflSlI, impact systems [iJI, etc. Analy zing them usually relies on the use of 
Lyapunov functions and input-to-state stability/nonlinear small-gain ideas [3y, ISOl [24l [23| , Lyapunov functionals |48l | 
or, more recently, another type of functionals, verifying certain boundary conditions, referred to as looped-functionals 
@, H Ell • When the impulses occur periodically, the system can be viewed as an LTI discrete-time system which 
can be studied using discrete-time Lyapunov theory. When impulses arrive at irregular times (aperiodic regime), 
the discrete-time system becomes time- varying and specific stability concepts should then be considered. The notion 
of dwell-time, i.e. the time between two successive discrete events defined as Tk := ifc+i — tfc, has been introduced 
early in the literature [471 134| and has been proven to be very useful for the analysis of switched systems. In the 
case of impulsive systems, dwell-times more specifically correspond to the times between two consecutive impulses. 
Impulsive systems can therefore be classified in terms of the properties of the considered sequence of impulse instants 
{tk}, and a relevant stability notion can be defined for that specific class of sequence of impulse times. When the 
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sequence of impulsive instants is arbitrary, i.e. T^ > 0, we talk about stability under arbitrary dwell-time, whereas 
stability under ranged dwell-time is defined for sequences verifying Tk G [TmimTmax]- Stability under minimum and 
maximum dwell-time address the cases Tk > T and Tk < T, respectively. 

Stability under dwell-time constraints can be analyzed in several different ways. Lyapunov approaches based on a 
separate worst-case convergence analysis (i.e. a-stability) of the distinct parts of the impulsive system [33, ll^ are 
very convenient to work with when dealing with uncertainties, or when control design is the main goal, principally due 
to their convexity properties. They have, however, been shown to be unable to characterize stability under ranged 
dwell-time, see e.g. Example 3 of Q or to yield very accurate estimates for dwcll-times |3ll |Tl[ since they may 
not capture the possible interplay between the continuous- and discrete-parts. Discrete-time approaches, however, 
exhibit much less conservatism, but are, in the present state-of-the-art, difficult to adapt to uncertain systems and 
aperiodic systems or to extend to control design, mainly due to a lack of convexity. Hybrid stability conditions 
(also referred to as "mixed stability conditions" in the following) consisting of coupled continuous-time and discrete- 
time criteria have been shown to yield more accurate estimates for minimum dwell-time for both linear switched 
systems [3l|,[21[ and linear impulsive systems [1, [I^. Considering such mixed conditions is, moreover, strengthened 
by the fact that the considered type of Lyapunov functions and type of conditions are simple instances of a general 
converse Lyapunov result proved in p38l | in the context of linear switched systems. The main difficulties in using 
hybrid conditions lie in the nonconvex dependence on the system matrices (due to the presence of a discrete-time 
condition), leading then to a difficult extension to both time- invariant and time- varying uncertain systems, and to 
control design. Looped- functionals [1, [l,[il[, which allow to convert discrete-time conditions into convex conditions, 
have made possible the extension of these mixed conditions to uncertain and aperiodic impulsive systems in a quite 
simple way, yet not computationally cheap. The problem of designing robust control laws with prescribed dwell-time 
using mixed stability conditions remains, however, a challenging and open problem. 

Periodic and aperiodic sampled-data systems, arising for instance in digital control jl7| or networked control systems 
[73 . [35I [7l| . are intimately connected to impulsive systems since any sampled-data system can be equivalently 
represented as an impulsive system. Several approaches have been developed to analyze sampled-data systems: 
discrete-time approaches |64. I28l 153 . 1371? ? ], input-delay a p proach es \6(k. 27 . 58tl2g . l43l .l4a|. robust analysis techniques 
[illlilli, impulsive/hybrid systems formulation [11 11 [11, llljpllini,!^ and the use of looped-functionals 
either considering directly the sampled-data system formulation [3, [S^] or the impulsive system formulation [^, |^, \1^ . 
As in the case of impulsive systems, the aforementioned approaches have both benefits and drawbacks. Indeed, while 
functionals or continuous-time robust techniques arc suitable for robust analysis thanks to the linear or convex 
dependence on the system matrices, they are, unfortunately, not very adapted to control design due to an increased 
number of decision matrices leading to cumbersome bilinear (nonconvex) terms in the design conditions; see e.g. 
[23, m, HH, d, [13 ■ Discrete-time approaches are, however, very well-suited for control design thanks to the presence 
of a single decision variable, i.e. the Lyapunov matrix, but may fail to yield tractable criteria when parametric 
uncertainties enter the picture, mainly due to the presence of exponential terms in the stability conditions; see e.g. 
[53l |. These exponential conditions are also limited to LTI sampled-data systems. 

The rationale for using mixed stability criteria for analyzing switched and impulsive systems [3l|, [2l|, [^ [13, [HI lies 
in the reduced (possibly vanishing) conservatism [68l | of the conditions, opposed to continuous-time results based on 
rates of convergence of Lyapunov functions, see e.g. [43, [13, [Sg . Hybrid stability criteria are therefore important to 
consider in order guarantee accuracy, but should be characterized in such a way that robustness analysis and control 
design remain possible. A first step towards such a result has been made very recently by using looped-functionals 
[a, [lOl [ill J ^ specific type of functionals that encode a discrete-time condition as a convex condition in the system 
matrices, a very suitable feature for robust stability analysis. However, the structure of the conditions prevents the 
derivation of tractable design criteria due to the presence of multiple decision matrices, inexorably leading to high 
computational cost and nonconvex terms in the synthesis conditions. The proposed approach, based on time- varying 
continuous-time discontinuous Lyapunov functions, combines features of the continuous-time and hybrid approaches, 
by leading to necessary and sufficient stability and stabilization conditions which arc convex in the system matrices 
and in the decision variables (Lyapunov and controller variables), together with a lower complexity than by using 
looped-functionals. 

The contribution of this paper is therefore manifold. First of all, a necessary and sufficient condition for stability of 
impulsive systems with periodic impulses are derived in Section 2 from the use of a specific discontinuous Lyapunov 
function. The advantage of the use of such Lyapunov functions lies in a reduced computational complexity over 
the use of looped-functional, while accuracy is mostly preserved. The periodic case is then extended to cope with 
aperiodicity in impulse arrival times (i.e. minimum, maximum, and ranged dwell-timcs) and time- varying parametric 
uncertainties. Necessary and sufficient results for discrete-time quadratic stability are provided, still with a reduced 



computational complexity. By relying on non-conservative algebraic manipulations, these results are, further, exactly 
adapted in Section 3 to quadratic (robust) stabilization using a particular class of state-feedback controllers. These 
stabilization results provide solutions to the impulsive system counterpart of the open problem of stabilizing linear 
switched systems with prescribed/optimized minimum dwcll-time; pUl penultimate paragraph, page 1916]. More 
concisely, quadratic stabilization with prescribed minimum, maximum or ranged dwell-times can be expressed as 
convex optimization problems. The approach is fully generic and can be applied to any linear impulsive system. 
Exploiting finally, in Section 4, the representability of sampled-data systems as impulsive systems, the results are then 
adapted to sampled-data systems. Convex necessary and sufficient conditions for quadratic stabilization of periodic 
and aperiodic uncertain sampled-data systems arc obtained. Examples and comparisons with several existing results 
are discussed in the related sections. 

Notations: The set of n x n (positive definite) symmetric matrices is denoted by (S"o) §". Given two symmetric 
matrices A, B, the inequality A >- (^)-B means that A — B ]s positive (semi) definite. Given a square real matrix A, 
the notation Sym(A) stands for the sum A + A^ . 

2 Stability analysis of periodic and aperiodic impulsive systems 

In the same vein as loopcd-functionals, the proposed approach is a hybrid approach yielding stability conditions 
convex in the system matrices. This allows us to take advantage of the structure of continuous-time criteria while 
benefiting of the accuracy of discrete-time criteria, as emphasized in past works such as [8J, [9|, |I0| . To develop further 
the approach, let us consider linear impulsive systems of the form 

x{t) ^ Ax{t), tj^tk 
x{t) = Jx^{t), t = tk 

where x € K" is the state of the system and x~{tk) stands for the left-limit of x{s) at s = tk, i.e. x~ (tk) = lim x{s). 

The system matrices A and J may be uncertain time- varying, this will be explicitly mentioned when this is the case. 
The sequence of impulse instants {ifcjfcgN, ifc > 0, is assumed to have positive increments Tk := tk+i — tk > e > 
that arc bounded away from 0. Defined as such, the sequence {tfejfeeN does not admit any accumulation point (we 
exclude then any Zeno motion) and grows unbounded. Note that the sequence of impulse instants may or may not 
depend on the state of the system. In the following, we will make no distinction between these two cases since impulse 
sequences will be solely characterized in terms of dwell-time constraints. 

We derive first necessary and sufficient conditions for stability of impulsive systems with periodic impulses which 
will serve as a groundwork for more advanced results on aperiodic impulsive systems with minimum, maximum and 
ranged dwell-times. These latter results are then finally extended to uncertain systems. Only the result on periodic 
impulsive systems is generalized to uncertain systems, all the others straightforwardly follow from the same lines. 

2.1 Impulsive systems with periodic im,pulses 

At first sight, the case of periodic impulsive systems, i.e. Tk = T, k G N, is of lesser interest since necessary 
and sufficient conditions for characterizing their stability already exist. The periodic case is, however, much more 
interesting than we may think, and this, for two reasons. First of all, this result allows us to present the main ideas 
in a simple and clear manner, easily generalizable to more elaborated cases such as time-varying impulsive systems. 
Second, the periodic result will turn out to be directly involved in aperiodic stability analysis under minimum and 
maximum dwell-times. 

Theorem 2.1 (Impulsive systems with periodic impulses) Let us consider the system (1) with periodic im- 
pulses, i.e. Tk ~ T, k IE N. Then, the following statements are equivalent: 

(a) The impulsive system (1) with T -periodic impulses is asymptotically stable. 

(b) The discrete-time transition matrix ^{T) :~ e J is Schu^_]. 



A matrix is Schur (or Schur stable) if all its eigenvalues lie in the unit disc. 



(c) There exists a matrix P G S"o such that the LMI 

J^e^'^Pe^^J -P<Q (2) 

holds or, equivalently, the quadratic form V{x) = x Px is a discrete-time Lyapunov function for the LTI 

discrete-time system Zk+i = e Jzk- 

(d) There exist a difjerentiahle matrix function R : [0, T] n- S", ^(0) >- 0, and a scalar e > Q such that the LMIs 

A^R{t)+R{t)A + R{t) ^0 (3) 

and 

J^R{0)J-R(f) + eI diO (4) 

hold for all r G [0, T] . 

(e) There exist a difjerentiahle matrix function S : [0,T] i— ;• §", S{T) >- 0, and a scalar £ > such that the LMIs 

A^S{t) + S{t)A - S{t) :< (5) 

and 

J^S{f)J-SiO) + eIdiO (6) 

hold for all t G [0, T] . 

Proof : Proof of (a) <^ (b) <^ (c): The equivalence between these statements is well-known and thus only sketched 
here. Since, the impulsive system is periodic then we can discretize the dynamics to get the equivalent discrete-time 
system x{tk+i) = e Jx{tk). The system (1) with T-periodic impulses is therefore asymptotically stable if and only 
if e J is Schur. This is also equivalent to the existence of a quadratic Lyapunov function V{x) = x Px that is 
pointwise decreasing along the trajectories of the discrete-time systems above. The proof is complete. 

Proof of (d) =^ (c): Assume (d) holds. Integrating (3) over [Q,T] implies that the LMIs 

e^'"'^R{f)e^^ - i?(0) d (7) 

and 

J^e^"^R{T)e^^J - J^R{0)J ^ (8) 

hold. From (4), wc have that R{f) y and J^i?(0)J ^ R{f) -el. Substituting then J^i?(0)J in (8) yields 

J^e''"^R{f)e^^J - R{f) ^-el (9) 

which therefore implies that (2) holds with P = R{T) >~ 0. The proof is complete. 

Proof of (c) => (d): The proof is structured as follows: first, we prove that (3) admits solutions regardless of the 
stability of the system, showing that this condition can be assumed to be satisfied without loss of generality. The 
second part of the proof consists of combining statement (c) with the solution set of (3) to prove that (4) holds. 

Assume (2) holds with P = R{T) >~ and some F ^ as 

J^e^'"^R{f)e^'^J - R{T) = -Y. (10) 

Since e J is Schur, then the above matrix equation admits the unique solution [3Q| | 

oo 

R{T) = Y, ( J^e-4"^) ' Y (e-4^ j) ' ^ 0. (11) 

i=0 



The set of all solutions R{t) to (3) can be defined as the set of solutions of the matrix equality 

A^R{t) + R{t)A + R{t) = -W{t), W{t) >r (12) 

where W{t) is a continuous function w.l.o.g. Given W{t), the unique solution to (12) is given by 

R{t) ^ e-^^"i?(0)e-^" - r e-^^^f"-^) VF(s)e-^("-"Ms, r G [0, f] (13) 

"'0 

where R{T) >- obeys (11). We have proved that (12) can be considered as fulfilled, independently of the stability 
of the system, which concludes the first part of the proof. The second part of the proof consists of deriving first, 
from expression (13), the equation 

e^^^R{T)e^^ ^ ^W{T) + R{0) (14) 

where W{f) = /J' e^'^^VF(s)e^Ms h 0. The above equality implies that i?(0) >- W{f) h since R{f) >- 0. 
Consequently, we have that R{t) >- for all r £ [0, T] since W{t) ^ is a nondecreasing function, i.e. W{t) ^ W{C) 
for any < T < C < f. 

Substituting, finally, the left-hand side of (14) in (10), we get that 

J'^R{0)J- R{f) = -Y + J^W(T)J. (15) 

Since W{s) and Y ^ are arbitrary, then we can choose them such that this then implies that —Y + J^W{T)J -< 
and thus that (4) holds. The proof is complete. 

Proof of (d) <^ (e): Assume (d) holds for some R{t), it is immediate to see that R{t) := S{T — t) solves (5) and 
(6). Reverting the argument proves the equivalence. <) 

The conditions stated in statement (d) can be understood as a non- increase condition, over each interval [tk,tk+i), 
of the time- varying discontinuous Lyapunov function 

Vd{x,T)=X^Q{T)x (16) 

where Q{tk + t) = Q{t), Q{t) G S", t G [0,Tfc), Q(0) >- 0, Q differentiable, and that verifies the boundary condition 

J^Q{0)J-Q-{f)+eI^O (17) 

where g"(r) = lim{g(s)}. 

A peculiarity of the proposed approach is that the matrices R{t) and S{t) do not need to be imposed to be positive 
definite over their domain of definition. Positivity over their domain is directly implied from the positivity of R{0) 
and S{T), and the LMI conditions in statements (d) and (e). These conditions, all together, indeed imply that R{T) 
is also positive definite, and thus that, by virtue of equation (13) that R{t) is positive definite on its domain. The 
case of S{t) is symmetric. 

There are several advantages of the conditions (3)-(4) of statement (d) (or conditions (5)-(6) of statement (e)) over 
condition (2) of statement (c). First of all, the conditions are convex in the system matrices A and J, allowing then 
for an immediate extension to the uncertain case. Further, the presence of a single decision matrix variable in the 
conditions tends to suggest the possibility of deriving tractable synthesis conditions. The compensation for these 
interesting convexity properties is the consideration of infinite dimensional feasibility problems, which may be very 
hard to solve. Several methods can be applied to make the feasibility problems finite-dimensional. The first one is to 
discretize the interval [0, T] and express the matrix R(t) as a piecewise linear function on each subintervals; see e.g. 
33, |l| . The second one relies on recent polynomial programming techniques [131 j such as sum of squares programming 
5J], which provide an efficient framework for solving such problems by restricting the matrix functions R{t) and 



S{t) to polynomial matrix functions. The package SOSTOOLS [53 together with the semidefinite programming 
solver SeDuMi [63l | supply the necessary material for solving such problems. The SOS program corresponding to 
statement (d) is expressed as follows 

Find polynomials R{t),M{t) G S" and e > 
such that M{t) is SOS 

-A^R{t) - R{t)A - R{t) - M{T)T{f - t) is SOS (18) 

i?(0) >- 

J'^R{Q)J -R{f)+el <0. 

The role of the SOS polynomial variable M is to incorporate the fact that the LMI (3) must be verified for all 
T £ [0,T] = {s £ R : s{T ~ s) > O}. This constraint incorporation can be justified using Putinar's Positivstellensatz 
|56l |. It is also very important to point out that the computation complexity is improved by the fact that R{t) docs 
not have to be specifically imposed to be positive definite over [0, T]. This constraint would indeed require the use of 
another SOS variable M2{t) in order to state the condition that R{t) — e I ~ M2{t)t{T — t) is SOS. This unnecessary 
complexity increase is circumvented by simply imposing i?(0) as positive definite. 

Still in a computational perspective, it seems necessary to compare the computational complexity of the conditions of 
Theorem 2.1 to the complexity of the looped- functional-based results of [l3| addressing the same problem. Assuming 
polynomial matrices R{t) £ §" (or S{t) G §") of degree dfi in Theorem 2.1 and a polynomial matrix Z{t) G S'^" of 
degree dz in [13, we have the following count of the number of variables: 

/, N /, -,n"-("- + 1) , ,>r /, N n(n + l) ,, _3n(3n+l) ,_,^, 

NcurrentidR) = {dR + 1)"^^ — - and Niooped{dz) = ^^ — - + {dz + 1)^ (19) 

for the current approach and the looped- functional approach of |10l |. respectively. We can immediately see that the 
number of variables for the looped- functional approach grows much faster with the system dimension n and the degree 
of the polynomial than with the proposed approach. It seems, however, important to stress that the expressions (19) 
should be understood as lower bounds on the actual computational complexity since additional variables are usually 
needed, e.g. to incorporate constraints. It will be illustrated in the examples that the proposed approach is able 
to obtain results that are very close to those obtained using looped-functionals with a much lower computational 
complexity, even ii d^ is usually larger than dz- A comparison will also be made with a discretization-based approach. 

2.2 Aperiodic impulsive systems 

Let us consider now that the system (1) is aperiodic, i.e. impulses arrive at irregular times. To this aim, we consider a 
ranged-dwell time constraint on the sequence of impulse instants, i.e. Tk G \Tmin, Tmax]- We then have the following 
generalization of Theorem 2.1: 

Theorem 2.2 (Ranged dwell-time) Let us consider the system (1) with a ranged dwell-time constraint, i.e. Tk G 
\T„iimTmax], fc G N. Then, the following statements are equivalent: 

(a) There exists a matrix P G S"o such that the LMI 

jTeA-ep^Aej_p^Q (20) 

holds for all 9 G [Tmin,Tmax]- 

(b) There exist a differentiable matrix function R : [0,Tmax] ^-^ §"; -R(O) >- 0, and a scalar e > such that the LMIs 

A^R{t)+R{t)A + R{t) ^0 (21) 

and 

J^R{Q)J - R{9) +el <Q (22) 

hold for all r G \^,Tmax\ and all 9 G [Tmin^Tmax]- 



Moreover, when one of the above statements holds, then the aperiodic impulsive system (1) with ranged dwell-time 
Tk G [Trmn,Tmax] «s asymptotically stable. 

Proof : The proof follows the same lines as the one of Theorem 2.1 

It is important to stress that, in the result above, the computational complexity of the second statement is much 
lower than if we had used conditions (3)-(4) of Theorem 2.1, statement (e). This follows from the fact that in 
statement (c), we would have required S{6) >- for all 6 G [0,Tmax], which is obviously much more complex than 
simply imposing -R(O) )^ in the present case. To pursue on the computational complexity analysis, we note that 
the conditions of Theorem 2.2 are more expensive than those of Theorem 2.1 due to the presence of the additional 
parameter £ [0,T,nax] in LMI (22). 

The next result concerns stability of impulsive system under minimum dwell-time, i.e. T^ > T for all k E N. This 
stability concept has been extensively studied in the past, see. e.g. [3^ d, [lO, 11] and references therein. 

Theorem 2.3 (Minimum Dwell-Time) Let us consider the system (1) with a minimum dwell-time constraint, 
i.e. Tfc > T, fc e N. Then, the following statements are equivalent: 



(a) There exists a matrix P G §"q such that the LMIs 

A^P + PA<Q (23) 

and 

J^e^^'Pe^'j-P^O (24) 

hold for all 9 G [T, oo). 

(b) There exists a matrix P e S"g such that the LMIs 

A^P + PA<Q (25) 

and 

jT^A-Tp^ATj_p^Q (26) 

hold. 

(c) There exist a differentiable matrix function R : [0, T] M- §", ^(0) >- 0, and a scalar e > such that the LMIs 

A^RiO) + R{0)A ^ (27) 

A'^R{t)+R{t)A + R{t) ^0 (28) 

and 

J'^R{0)J - R{f) + e / r< (29) 

hold for all r G [0, T] . 

(d) There exist a differentiable matrix function S : [0, T] M- §", S{T) >- 0, and a scalar £ > such that the LMIs 

A^S(T) + S{f)A -< (30) 

A^S{T) + S{T)A~S{T)diO (31) 

and 

j'^S{f)J-S{0) + eI :<0 (32) 



hold for all t G [0, T] . 

Moreover, when one of the above statements holds, the impulsive system (1) is asymptotically stable under minimum 
dwell-time T, i.e. for any sequence {tk}keN such that T^ > T. 

Proof : Proof of (b) <^ (c) <^ (d): It follows from Theorem 2.1. 

Proof of (a) ^ (6): Immediate. 



Proof of (b) ^ (a): This has been proved in [^, \T^, this is recalled for completeness. The key idea is to prove 
that (25) and (26) together imply that (24) holds for aU t € [f, cx)). Let €{9) := e^^^Pe-^^ then (24) holds for all 
e [f, oo) if and only if j'^ C{e)J -P <Q holds for aU 9 £ [t, cx)). A Taylor expansion around 9 = 9o>f yields 



£(610 + S)- C{9o) = Se^^^°[A^P + PA]t 



Ae„ 



o{S). 



(33) 



Since (25) holds, then the right-hand side is negative definite for all 9q > T and all 5 > 0, thus we have that 
C(T + e) ^ C{f) for all e > 0. Therefore, j'^ C{f)J - P ^ implies that j'^ C{e)J - P ^ holds for all 9 > T, 
which concludes the proof 

It is important tot note that the above theorem straightforwardly extends to time- varying systems depending ex- 
plicitly on time and/or time- varying parameters by simply using the fundamental-solution and the state-transition 
matrices instead of exponentials. The variational argument to prove the equivalence between statements (a) and (b) 
remains also valid. 

Remark 1 Similarly to as in Q, Hul l, a maximum dwell-tim,e result can be obtained by simply reverting the inequality 
sign in LMIs (23) , (25), (27) and (30). In such a case, the concluding statement changes to: "The aperiodic impulsive 
system (1) is asymptotically stable under maximum dwell-time T , i.e. for any sequence {tk\k&i such that Tk S [e, T] 
for any e > . " 

2.3 Examples 

The conditions stated in Theorems 2.1, 2.2 and 2.3 are infinite-dimensional feasibility problems. In order to enforce 
them efficiently, the sum-of-squares programming package SOSTOOLS [55] and the semidefinite programming solver 
ScDuMi [631 ^-^"^ used. Suitable matrix functions R or S such that the conditions of Theorems 2.1, 2.2 and 2.3 are 
feasible are then searched within the set of matrix polynomials of fixed (and chosen) degree, du say. In this case, 
the matrix function R{t) is chosen as R{t) = X]i=o ^i'^^^ ^i ^ ^'N ^^^i ^^ this regard, its derivative is simply given 
by the polynomial R{t) = J2i=i iRiT^^^ which can be easily inserted in the SOS conditions. In the examples below, 
the number of variables is identified as the number of variables declared by SOSTOOLS when defining the matrix 
decision variables, i.e. the Lyapunov matrix P{t) and the SOS variables Mi(r)'s for constraints incorporation. 

Example 1 (Ranged dwell-time) Let us consider the system (1) with matrices 0, \luJ 



A = 



"-1 O.l' 


, J ^ 


"l.2 


o" 


1.2 







0.5_ 



(34) 



By computing the eigenvalues of e J, this system can be easily shown to be stable with T -periodic impulses whenever 
T £ [0.1824,0.5776]. Stability of this system in the aperiodic case cannot be directly addressed using methods based 
on rates of convergence of Lyapunov functions since both matrices are unstable, see e.g. fSa . [3/. Note, however, that 
this system can be analyzed using methods from f2ji] . Using the ranged dwell-time stability conditions (21)-(22) of 
Theorem 2.2, the same bounds are retrieved with a matrix polynomial R of order 6, showing then tightness of the 
obtained numerical values in the aperiodic case. For comparison purposes, the same numerical result is obtained 
in flQ ] using a looped- functional of degree dz = 3. SOSTOOLS, however, declares 14-9 variables for the current 
approach, whereas for looped-functionals 2806 variables are involved. The execution time, on an Intel il, 2.70 Ghz, 
4GB of RAM, including pre- and post-processing is about 1 second whereas it is approximately of 15 seconds for the 
looped functional approach of WXj . 

Example 2 (Minimum dwell-time) Let us consider the system (1) with matrices 0, [73 / 



A = 



"-1 


o" 


, J = 


"2 1' 


1 


-2 




1 3 



(35) 



Since A is Hurwitz, the minimum dwell-time result stated in Theorem 2.3 can be applied. Using conditions (25)- 
(26), we get the minimum dwell-time T = 1.1405. T/ie same value for the minimum dwell-time is obtained using 





order of R 


-/ min 


-i max 


Theorem 2.2, (b) 


2 
4 
6 


0.1834 
0.1824 
0.1824 


0.4998 
0.5768 
0.5776 


Periodic case 


- 


0.1824 


0.5776 



Table 1 

Estimates of the admissible range of dwell-times for the 

aperiodic system of Example 1 





order of R 


T ■ 

-*- mm 


Theorem 2.3, (c) 


2 
4 
6 


1.1883 
1.1408 
1.1406 


Theorem 2.3, (b) 


- 


1.1406 


Periodic case 


- 


1.1406 



Table 2 

Estimates of the minimum dwell-time for Example 2 



conditions (27)-(28)-(29) with a polynomial matrix R of order 6; see Table 2. Using the looped-functional approach 
of flQ l. this numerical result is obtained by using polynomials of degree dz ^ 3 (i.e. ^12 variables), whereas the 
current approach involving polynomials of order 6 only requires 85 variables; approximately a five-fold decrease. 
The execution time, on an Intel i7, 2.70 Ghz, ^GB of RAM, including pre- and post-processing is about 0.5 second 
whereas it is approximately of 1.5 seconds for the looped functional approach offlQL For comparison, we also consider 
a discretization scheme Q/ where R{t) is expressed as a piecewise linear function on [0, T] which is subdivided in N 
subintervals. For a fair comparison, we select N = 28, which gives 87 variables, and we get 1.1919 as the computed 
bound on the minimum dwell-time. The computation time is 1.2 seconds. Thus we can see that, on this example, the 
SOS approach perform better with a comparable number of variables. Note, moreover, that the number of constraints 
involved in the discretization is larger than the one considered in the SOS program as well. 

Example 3 (Maximum dwell-time) Let us consider the system (1) with matrices 0, [7fl / 



A 



\ ^ ^1 






, J ^ 


[-1 2J 





1/2 
1/2 



(36) 



This system is a candidate for the maximum dwell-time result discussed in Remark 1. Conditions (25)-(26), where the 
inequality sign of (25) is reversed, yield the maximum dwell-time T = 0.4620. Using then conditions (27)-(28)-(29) 
with a sign reversal on (27), we obtain the maximum dwell-time values summarized in Table 3. The value 0.4620 is 
retrieved using a matrix polynomial R of order 6. Once again the computational complexity of the approach is much 
lower than by using looped- functionals; see Example 2. 





order of R 


-i max 


Theorem 2.3, (c) 


2 
4 
6 


0.3999 
0.4615 
0.4620 


Theorem 2.3, (b) 


- 


0.4620 


Periodic case 


- 


0.4620 



Table 3 

Estimates of the maximum dwell-time for Example 3 using The- 
orem 2.3 and Remark 1. 



2.4 Robustness results 



order of R 


-t max 


2 

4 


0.1028 
0.1071 


- 


0.1072 (gridding) 



Table 4 

Estimates of the maximum dwell-time for Example 4 using the 
robustification of Theorem 2.3 (using Theorem 2.4) and Re- 
mark 1. 



All the previous results can be robustified to account for parametric uncertainties affecting A and J. To this aim, 
let us consider now that the matrices of the system (1) are uncertain, possibly time- varying, and belonging to the 
following polytopes 

A(E A:=co{Ai,...,An} and J e J' := co{Ji,..., Jjv} (37) 

where co{-} is the convex- hull operator. Before stating the main results, it is necessary to introduce the state- 
transition matrix $(•), which corresponds to system (l)-(37), as 



d$(s) 
ds 



= ff^A,(s)Aj$(s), $(0)=/ 



(38) 



where A(s) G Ajv := {S, G K>q : ||^||i = l} is sufficiently regular so that solutions to (38) are well-defined, e.g. in a 
Carathcodory sense. Associated to this transition matrix, we define the set $y as 

^f := {$(T) : $(s) solves (38), A(s) 6 Ajv, se [0,f]} . (39) 

This set corresponds of all possible transition matrices $(r) obtained for all possible trajectories of the uncertain 
parameters A. Note that the set $y is strongly nonconvex and is difficult to compute exactly. This intricate structure 
illustrates the inherent difficulty in considering uncertain systems in a discrete-time setting. By reformulating the 
discrete-time conditions in terms of conditions (3) and (4), this difficulty is circumvented and discrete-time stability 
results can be efficiently robustified. For conciseness, only the robustification of Theorem 2.1 will be discussed. 
Robust versions of Theorems 2.2 and 2.3 can be obtained in the same way. 

Theorem 2.4 (Periodic impulses) Let us consider the uncertain (time-varying) impulsive system (l)-(37) with 
T -periodic impulses, i.e. T^ ^ T, fc G N. Then, the following statements are equivalent: 

(a) The uncertain (time-varying) impulsive system (l)-(37) with T-periodic impulses is quadratically stabla^. 

(b) There exists a matrix P G S"o such that the LMI 

J^^iJ'^pq)J - P ^0 (40) 

holds for all (^, J) G #y x J^. Equivalently, the quadratic form V{x) = x Px is a discrete-time Lyapunov 
function for the uncertain time-varying discrete-time system Zk+i = ^kJzk, for all (^fe, J) G ^j- x J^ . 

(c) There exist a differentiable matrix function R : [0, T] l-^' §", ^(0) >- 0, and a scalar e > Q such that the LMIs 

AJR{t) + R{t)A, + R{t) r< (41) 

and 

jfi?(0)J, -i?(T)+e/r<0 (42) 

hold for all r G [0, T] and all i — 1, ... , N . 

(d) There exist a differentiable matrix function S : [0,r] i— 5- §", S{T) >- 0, and a scalar £ > such that the LMIs 

AJS{t) + SiT)A, - S{t) ^ (43) 

and 

J^SiT)J,~SiO) + eI^O (44) 

hold for all r G [0, T] and all i — 1, ... , N . 

N N 

Proof : Let us introduce the matrices ^(A) = 2, ^lAi and J(C) — 2, QJi^ ^^C ^ ^n- The proof then simply relies 
on convexity arguments. Indeed, the LMI 

A{X{T)fR{T) + R{t)A{X{t)) + R{t) ^ (45) 

holds for all A(r) G An and aU r G [0, f] if and only if 

AfR{T) + R{t)A, + R{t) < 
holds for all r G [0, f] and alH = 1, . . . , A^. Similarly, using the fact that i?(0) >- 0, the LMI 

J{CfR{Q)J{C) - R{T) + e / ^ 0, C e Ajv (46) 



Quadratic stability of a linear uncertain system is defined here through the existence of a common quadratic Lyapunov 
function (i.e. independent of A in the present case) for the uncertain system; see e.g. [40]. 
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is convex in C, and therefore holds for ah C G A^v if and only if 

jfi?(0)J, -i?(f') + e/r<0 

holds for all i = 1, . . . ,N. The rest of the proof follows from the equivalence between (45)- (46) and (40) using the 
same arguments as for proving Theorem 2.1. Note that the proof of Theorem 2.1 readily extends to the case of time- 
varying matrices where exponential terms simply have to be replaced by state-transition matrices belonging to ^f. <0> 

Example 4 Let us consider system (l)-(37) with the polytopes flQ] 



A 



' 1 3 




"3 l' 


L j = 


= CO < 


"l/4 


' 




"2/3 





-1 2 




-2 4 


1 


. 


_ 


2/3 




_ 


1/2 



(47) 



This system is the uncertain version of the system treated in Example 3. By gridding the conditions (25) -(26) where 
the first one has inequality sign reversed, we obtain the optimistic value (upper-hound) 0.1072 on the robust maximum 
dwell-time. Applying, on the other hand, the robust version of the maximum dwell-time result^, we obtain the results 
of Table 4 where we can see that the value 0.1072 is almost recovered whenever R is chosen as a matrix polynomial of 
order 4. For comparison, the same result is obtained with a looped-functional of degree dz = 2 in flOl L However, the 
current approach involves 87 variables, whereas 335 variables are declared by SOSTOOLS for solving the conditions 
obtained using looped- functionals. The execution time, on an Intel i7, 2.70 Ghz, 4GB of RAM, including pre- and 
post-processing is about 1 second whereas it is approximately of 4 seconds for the looped functional approach of yWj . 

3 Stabilization of periodic and aperiodic impulsive systems 

It is now shown that the current framework can be efficiently and accurately used for control design, unlike loopcd- 
functionals and several approaches based on Lyapunov functionals. To this aim, let us consider the impulsive system 



x{t) = Ax{t)+B,u^{t), t^tk 
x(t) = Jx-{t) + BdUd{t), t^tk 



(48) 



where Uc € W^" and Ud G K™'' are the control inputs. The following class of state-feedback control laws is considered: 



Uc{tk + r) = K,{T)x{tk + t), re [0,Tfe), A; G N 
■Ud{tk) = Kdx'^itk) 



(49) 



where the continuous control law is time-varying and the discrete one is time-invariant. The purpose of this section 
is therefore to provide tractable conditions for finding suitable Kc : [0,T) i-^ ^m^xn ^^^^ ^^ ^ j^m^xn ^^^j^ ^j^^^ ^j^^ 
closed-loop system (48)-(49) is asymptotically stable. 

In the following, exact conditions for nominal stabilization of impulsive systems with periodic impulses and minimum 
dwell-time arc provided. The non-conscrvativeness of the conditions is strongly related to the considered control law: 
choosing, for instance, a constant gain for the continuous control law would have resulted in much less tractable 
conditions, yet possible to handle conservatively using various relaxation techniques. As for stability analysis re- 
sults, stabilization results on ranged and maximum dwell-times readily follow; they are thus omitted. Extensions to 
uncertain systems are immediate and omitted as well. 

3.1 Periodic impulses case 

The following result gives constructive conditions for designing a control law of the form (49) for impulsive systems 
with T-periodic impulses, i.e. T^ = T. Suitable controller gains can, indeed, be directly extracted from the solutions 
of the sum-of-squares feasibility problems stated in the following result: 



Not explicitly derived in this paper but consisting of the robustification of Theorem 2.3, and accounting for Remark 1. 
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Theorem 3.1 (Impulsive systems with periodic impulses) The following statements are equivalent: 



(a) There exists a control law (49) such that the impulsive system (4-8)-{A9) with T-periodic impulses is asymptoti- 
cally stable. 

(b) There exist a differ entiable matrix function R : [0,T] i— > §", RiT) )^ 0, a matrix function Uc '■ [0,T] i— > R™<:^"^ 
a matrix Ud G K"''*'^" and a scalar e > such that the LMIs 



Sym[Ai?(r) + B,Uc{t)] - R[t) < 



and 



-<0 



-R{Q) + eI JR{T) + BdUd 

• -R{f) 

hold for all r G [0, T] . In such a case, suitable matrices for the control-law (49) are given by the expressions 



(50) 
(51) 



Kc{t) = Uc{t)R{t)-' 
Kd = UdR{f)-\ 



(52) 



(c) There exist a differentiable matrix function S : [0,T] i-> §", S{0) >- 0, a matrix function Uc : [0,T] i-^ jg'ncXn^ 
matrix Ud G R"'''^" and a scalar e > such that the LMIs 



Sym[AS{T) + S,C/,(t)] + S{t) ^ 



and 



-<0 



-S{T) + eI JS{0) + BdUd 

• -S{Q) 

hold for all r G [0, T] . In such a case, suitable matrices for the control law (49) are given by the expressions 



(53) 
(54) 



Kcir) = [/e(r)^(T)-i 
Kd = UdS{0)"\ 



(55) 



Proof : What has to be proven is the exactness of the stabilization conditions (53)-(54). By performing a congruence 
transformation on (53) with respect to S := S~ we get that 



Sym 



SiT)iA + BcKciT)) -SiT)^0 



(56) 



where we used the facts that Kc{t) ~ Uc{t)S{t) and S{t)S{t)S{t) = —S{t). Looking now at the LMI (53), we can 
easily sec that it is equivalent to 



(J + BdKdf 



-S{T) 
S{0) 



{J + BdKdY 



^0. 



(57) 



Noting then that the central matrix has n positive and n negative eigenvalues, and that the outer-factors arc of rank 
n, then the dualization lemma [H, [13 applies, and we get the equivalent LMI 



(J + BdKdfS{f){J + BdKd) - S{0) < 



(58) 



and thus that 

(J + BdKdfSiT){J + BdKd) - S{0) +e'I^O (59) 

for some e' > 0. Noting finally that the conditions (56)-(59) are identical to (5)-(6) proves the result. The equivalence 
follows from the losslessncss of the manipulations. Statement (b) can be proved along the same lines. <0 
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Remark 2 Whenever stabilization with ranged dwell-time is considered, the conditions of statement (h) should he 
considered in order to get a controller matrix K^ that is independent of Tk since the associated decision variable 
Ud is coupled to Kd via the matrix S{0). Statement (c) will, on the other hand, lead to a controller matrix Kd that 
depends on Tk since Ud and Kd are coupled via S(Tk). 

3.2 Stabilization under minimum dwell-time 

The stabilization under minimum dwell-time is slightly more complicated since the controller gain K(.{t) in (49) 
must remain bounded as r — >■ cxd. In the best case, it should converge to a finite value. A way for solving this difficulty 
is to consider the following controller gain 

\/i,(T) ifT€ [T.Tt) 

where T^ > T, fc G N and Kc{t) is some matrix function to be determined. This specific structure for the control 
law, as it will be emphasized later, arises naturally from the structure of the minimum dwell-time stability conditions 
and will be shown to be non-restrictive. Again the matrices of the controller can be extracted from the solutions of 
the feasibility problem stated in the following result: 

Theorem 3.2 (Minimum dwell-time) The following statements are equivalent: 

(a) There exist matrices P G §"o; ^^d £ R™''^" and a matrix function Kc '■ [^,T] h^ R^cXn ^^^^ ^^^^ ^^g ^jiatrix 
inequalities 

{A + BJ<,{T)fP + P{A + BJ<,{T)) -< (61) 

and 

(J + BdKdf<i>{ffPHf){J + BdKd) -P<Q (62) 

holds where $ : [0, oo) M> M'"^" is the transition matrix defined as 

^$(r) = [A + B,K,{T)Mr), r > ^^^^ 

$(0) = /. 

(b) There exist a differentiable matrix function S : [0,r] M- §", S{T) ;^ 0, a matrix function Uc ■ [0,T] n-R™''^", 
a matrix Ud G R'"<i^" and a scalar e > Q such that the LMIs 

Sym[ASiT) + B,Uc{f)] < 0, (64) 

Sym[^S'(r) + B^Udr)] + S{t) < (65) 



and 

^ -S{Q) + eI JS{T) + BdUd 



^ (66) 



-S{T) 
hold for all r £ [0,r]. In this case, suitable controller gains are retrieved using 

Kc{r) = U,{T)SiT)~' 

- 1 (67) 

Kd = UdSiT)-\ 

Moreover, in such a case, the closed-loop system (48) -(49) -(60) is asymptotically stable with minimum dwell-time T. 

Proof : The first thing that has to be proven is the fact that statement (a) implies that the closed-loop system is 
stable with minimum dwell-time T. The equivalence between (a) and (b) follows from Theorem 3.1 and the changes 
of variables Uc{t) — Kc{t)S{t) and Ud = KdS{T). Let us prove then that statement (a) implies that the closed-loop 
system is stable with minimum dwell-time T. Two possible scenarios: 1) cither impulses always arrive in finite-time, 
i.e. T < Tk < oo; or 2) impulses stop at some point. 
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Case 1): In the first case, we need to show that the Lyapunov function V{x) = x^Px evaluated at times tk and 
along the trajectories of the closed-loop system (48)-(49)-(60) is pointwise decreasing. This is the case whenever the 
LMI 

(J + BaKaf<^{efP<^{e){J + B^Kd) ~P<Q (68) 

holds for all finite 9 e [f,oo). By using the fact that $(f + 5) = e(^+'^<=-^'=(^»'^$(T) for all (5 > and the same 
argurnents as in the proof of Theorem 2.3, we have that conditions (61) and (62) implies that (68) holds for all finite 
9 € [T, do). The convergence of the trajectories of the impulsive system to simply follows from the boundedness 
and continuity of the function V{x{t)) on every interval {tk,tk+i). This proves case 1). 



Case 2): When impulses stop after some time, tk say, then only the continuous regime is active and the system 
becomes a standard continuous-time system whose asymptotic stability is governed by the long-term behavior of the 
continuous-time matrix A -f BcKc{t). Since we have 



A + B,K,{t) 



A + BcKc{t) ifrG[0,r) 
A + BcKciT) iiT>T 



(69) 



then, by virtue of the condition (61), the Lyapunov function V{x) = x Px starts decreasing at t^ + T, the latest, 
and continuously decreases to 0, implying then asymptotic stability of the system. This proves case 2). {> 

In the above result, we can clearly see that the structure of the control law fits exactly the structure of the conditions, 
which allows us to obtain lossless results. Without this maintain of the value of K{t) to K{T) for all r G \T,Tk), 
deriving a minimum dwell-time stabilization result would have been much trickier. As a concluding remark on 
stabilization under minimum dwell-time, we note that the proposed control law just needs to be computed offline 
and is easy to implement. 

Example 5 Let us consider the system (48) with matrices 



A 



"i o' 


, B = 


"l" 


1 2 








and J = 



1 1 
1 3 



(70) 



Note that this system has both unstable flow and jumps. If we therefore assume that Kd = 0, then the system cannot 
be stabilized for arbitrary dwell-times since the jumps are destabilizing. Thus, we want to compute Kc{t) such that 
the minimum dwell-time is, at most, T = 0.1. Invoking Theorem 3.2, statement (b), with polynomial matrices Uc{t) 
and P{t) of order 1, we find the controller 



Kc{t) 



d{r) 



1.4750481 + 3.2714889r - 41.011914r" 
3.9063911 - 1.6733059r - 37.472443r2 



(71) 



where d{T) — —0.19767438 + 0.78454217r 4- 7.6562219r . State-trajectories of the closed-loop system, for some 
randomly generated impulse-times satisfying the dwell-time constraint, are depicted in Fig. 1. We can clearly see 
that the controller stabilizes the system. In terms of computational complexity, only 27 variables are defined by 
SOSTOOLS. 

4 Application to aperiodic uncertain sampled-data systems 

It is well-known that sampled-data systems can be reformulated as impulsive systems. On the basis of this refor- 
mulation, all the results developed in the previous sections therefore apply. It is hence possible to obtain robust 
stabilization results for aperiodic uncertain sampled-data systems with inter-sampling times in a range, providing 
then a solution to this challenging problem. In terms of robust stability analysis, the proposed approach is computa- 
tionally less expensive than those based on looped- functionals; see [13, |6^- Input-delay based methods [23, EJ are 
less expensive but less accurate as well, and cannot be efficiently used for control without drastic simplifications. 
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2.5 



Fig. 1. State-trajectories of the closed- loop system (70). 

Discrete-time approaches perform very well for stabilization, see e.g. |53l |. but are, until now, difficult to extend to 
uncertain systems due to the presence of exponential terms. The advantages of the proposed approach lie in the 
derivation of accurate convex robust stability and stabilization results together with a lower complexity than by 
using looped-functionals [10, Hol • 

4-.1 Preliminaries 

Let us consider here the continuous-time system 

x{t) = Ax{t) + Buit) (72) 

where x G K" and u e M™ are the state of the system and the control input, respectively. The control input is 
assumed to be computed from a sampled-data state-feedback control law given by 



L{t) = Kix{tk)+K2u(tk^l), t e [tfc,tfc+i) 



(73) 



where Ki € K^x" and K2 G K™^™ arc the control gains to be determined. Note that the control-law, viewed as 
a discrete-map from x to u, is BIBO-stable if and only if K2 is Schur. Above, the sequence of sampling instants 
{^fcjfceN is assumed to be strictly increasing and unbounded, i.e. t^ -^ 00, excluding therefore any Zeno behavior. 

The sampled-data system (72)- (73) can be equivalently reformulated as the following impulsive system 



'x{t) 


= 


'a b 






'x{t) 
_z{t)_ 



t^tk 



'x{t) 




' I 


0' 




'x-{t) 


[z{t)\ 




Ki 


A'2. 




[^-m 



(74) 



, t^tk 



where z G M™ is an additional state containing the value of the held control input at any time, i.e. z{t) — u{tk), 
t G [tfc,^/£+i)- For convenience, we also decompose J as J = Jq + BqK where 



Jn = 



7 0' 


, Bo = 


"0" 







/ 



and K = 



K, Ko 



(75) 
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4-2 Stabilization of periodic and aperiodic sampled-data systems 

Stabilization of sampled-data systems being the most interesting problem, we will therefore focus on the stabilization 
of both periodic and aperiodic sampled-data systems. As in previous stabilization results, a suitable stabilizing 
controller gain K can be extracted from the solutions of a feasibility problem. 

Theorem 4.1 (Periodic sampled-data systems) The following statements are equivalent: 

(a) There exists a gain K G ]]jmx(m+n) ^^^^ ^^q,^ ^^g T-periodic sampled-data system (72)-(73) is asymptotically 
stable. 

(b) There exist matrices P G §"^" and U G R™x(»+™) such that the LMI 



-P e^^iJoP + BoU) 
* -P 



-<0 



(76) 



holds. Moreover, in such a case, a suitable stabilizing controller gain is retrieved using K ~ UP 



(c) There exist a differentiable matrix function R : [0,T] i— S- 



:^n+m 



> such that the conditions 
and 



S{0) yO, a matrix Y G M'»x(»+™) and a scalar 



A{t)S{t) + S{T)A{Tf + S{t) :< 



-<0 



-S{T)+sI JoSiO)+BoY 
* -S{Q) 

hold for all t G [0, T] . Moreover, in such a case, a suitable stabilizing control gain is retrieved using K = Y S{Q)~ . 



(77) 
(78) 



Proof : The proof relies follows the same lines as the ones of the proofs of the stabilization results developed for 
impulsive systems. <^ 

Remark 3 Interestingly, it is also possible to impose K2 ~ without introducing any conservatism. Such a controller 
can indeed be designed by simply imposing the value to n x m right-upper block of the matrix S{0). The reason 
why this equality constraint is non-restrictive lies in the fact that e J is a block lower triangular matrix, and that 
it is well-known that stability of cascade systems (that are represented in terms of a block triangular matrix) can be 
exactly characterized by block diagonal Lyapunov functions. 

When the system is aperiodic, the following result can be obtained: 

Theorem 4.2 (Aperiodic sampled-data systems) The following statements are equivalent: 

(a) There exists a control law of the form (73) that quadratically stabilizes the system (72) for any aperiodic sampling 
instant sequence {tk} such that Tk G [TmimTmax]- 

(b) There exist matrices P G §"^™ and U G R"x("+™) such that the LMI 



-P e^^JoP + SoC/) 
* -P 



-<0 



(79) 



holds for all 6 G [TmimTmax]- Moreover, when this statement holds, a suitable stabilizing control gain can be 
obtained using the expression K = UP~ . 
(c) There exist a differentiable matrix function R : [Q,Tmax\ ^ S"+"\ S{0) y 0, a matrix Y G M™x("+™) and a 
scalar e > such that the conditions 



AiT)SiT) + SiT)A{Tf + Sir) ^ 



(80) 
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and 

' -S{9)+eI Jo + BnY 

• -S{0) 

hold for all r G [0,T.max] o,nd all G \T,n.in,Tmax\- Moreover, when this statement holds, a suitable stabilizing 



<0 



(81) 



control gain can be obtained using the expression K = YS{0) 
Example 6 Let us consider the sampled-data system (72) with matrices 



A 



1 

-0.1 



and B 




0.1 



(82) 



Fixed control law: Assume first that the control law is given as in 0/ by Ki = —3.75 —11.5 o-nd K2 = 0. Results 

in the aperiodic case (Tmin has been set to 0.001) are summarized in Table 5 together with some comparisons with 
previous ones based on Junctionals. The proposed approach yields results that are very close to the looped-functional 
approach developed in ffwl l together with a reduced computational complexity. The semidefimte program generated by 
SOSTOOLS involves, when R or S is of degree 4, 192 variables whereas the approach of fOul involves 1256 variables 
when using a polynomial of order 3. The execution time, on an Intel il, 2.70 Ghz, 4GB of RAM, including pre- and 
pos t-processing is about 1 second whereas it is approximately of 4-46 seconds for the looped functional approach of 

Control design: Assume now that the control gains Ki and A'2 have to be determined such that the closed-loop 
system is stable for any inter-sampling times in [TmimTmax]- Applying then Theorem 4-2, we obtain the results 
gathered in Table 6 where, following Remark 3, K2 = has been imposed in the three last scenarios. We can see that 
the computed controllers involve numerical values with reasonable magnitude. The involved number of variables for 
degrees of R or S equal to 2 and 3 are 201 and 432, respectively (one variable less when A'2 = Q). 

Example 7 Let us consider the sampled-data system (72) with matrices ]6_ 



A = 








-0.9 



, B = h, Ki = 



-1 



and K2 = 0. 



(83) 



The results in the aperiodic case are summarized in Table 5. We can see that results using looped-functionals from 
wU l can be almost recovered when R or S is of degree 4, which corresponds to 338 variables. The looped-functional 
approach of wOl l considers a matrix function of order 5, which corresponds to 3414 variables. The execution time, 
on an Intel il, 2.10 Ghz, 4GB of RAM, including pre- and post-processing is about I.I4 seconds whereas it is 
approximately of 15.48 seconds when using looped-functionals wOi l. 

Example 8 Let us consider the sampled-data system (72) with matrices 0/ 



A 



1 
-2 0.1 



and B 



(84) 



borrowed from the time-delay system literature fSs l. Assuming the control law Ki =10 and K2 = 0, we get 

the results of Table 5. We obtain a result very close to the one of fOul l using a matrix R or S of degree 6, which 
corresponds to 330 variables. The result of JOu l considering a matrix polynomial of order 5 involves 3414 variables. 
The execution time, on an Intel il, 2.70 Ghz, 4GB of RAM, including pre- and post-processing is about 1.12 seconds 
whereas it is approximately of 15.34 seconds for the looped functional approach of fOu]. 

4.3 Robust stabilization of periodic and aperiodic sampled-data systems 

Results on robust stabilization of sampled-data systems are straightforward extensions of Theorems 4.1 and 4.2; 
they are omitted for brevity. Only the following example is discussed: 
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Table 5 

Estimates on the minimum and/or maximum sampling period for the systems of Examples 6, 7 and 8 - Aperiodic case. 





order of 7? or S 


System (82) 


System (83) 


System (84) 


J- max 


J- min 


J- min 


J- max 


Th. 2.1 


2 
4 
6 


1.0851 
1.7279 
1.7252 


2.1496 
3.2480 
3.1933 


0.4 
0.4 
0.4 


0.8055 
1.6316 
1.8270 


[27] 
[48] 
[26] 
[43] 
[59] 


- 


0.869 
1.113 
1.695 
1.695 
1.723 


0.99 
1.88 
2.03 
2.53 
2.62 


0.400 


1.251 


[60] 


3 

5 


1.7294 
1.7294 


3.218 
3.269 


0.4 
0.4 


1.820 
1.828 



Table 6 

Control design results for system (82) using Theorem 4.2 



J- min 


J- max 


Ki 


A'2 


order of R 


0.001 


6 
10 
50 




-0.2185 -1.3276 
-0.1145 -0.8088 
-0.0202 -0.1560 




0.0022 
-0.0024 
-0.0030 


2 
2 
2 


0.001 


10 
20 
50 




-0.0310 -0.3222 
-0.0074 -0.0811 
-0.0259 -0.2726 









3 
3 

4 



Table 7 

Control design results for system (85) using Theorem 4.2 



s 


J- min 


J- max 


A'l 


K2 


order of R 


5 
5 


0.001 


10 
20 




-0.0757 -0.7306 
-0.0411 -0.3835 




-0.0006 
-0.0022 


2 
2 


20 
20 


0.001 


10 
20 




-0.0578 -0.5560 
-0.0339 -0.3121 




-0.0025 
-0.0019 


2 
2 



Example 9 Let us consider the uncertain sampled-data system (72) with matrices 



Ae A = CO- 



"0 


1 


J 


"0 


1 





-0.1 







-0.1_ 



and B 



(85) 



where S is a positive parameter. We then apply Theorem ^.2 to design robust state-feedback controllers for different 
values for S > and Tmax > 0. The results are summ,arized in Table 7 where we can see that the system can be 
stabilized for a quite wide range of values for the parameter 5 and the maximal sampling period T^ax ■ For a polynomial 
or order 2, the semideflnite program involve 231 variables. The execution time including pre- and post-processing is 
about 2.23 seconds on an Intel il, 2.70 Ghz, 4GB of RAM. 

5 Conclusion 

Discontinuous Lyapunov functions have been shown to be a powerful tool for analyzing stability of linear impulsive 
systems. The flexibility of the framework allows for easy extensions to aperiodic systems and uncertain systems 



18 



together with a reduced computational complexity with respect to looped-functional based approaches. The presence 
of a unique decision matrix permits the derivation of convex synthesis conditions without any simplifying nor 
restrictive assumptions. This method has been shown to provide a solution to the open problem of deriving convex 
conditions for robust stabilization of aperiodic uncertain linear impulsive with minimum, maximum and ranged 
dwell-times. 

As a byproduct, the approach can also be applied to sampled-data systems by exploiting the generic representability 
of sampled-data systems as impulsive systems. By transitivity, the approach has been shown to provide a solution to 
the open problem of deriving convex conditions for robust stabilization of aperiodic uncertain time- varying sampled- 
data systems. 

Addressing the cases of homogeneous polynomial Lyapunov functions pOl [2ll | and nonlinear impulsive systems 
are possible generalizations of the approach. Notably, homogeneous Lyapunov function can be used to vanish the 
conservatism of the approach for linear systems, following to [H, [2l|. There is also hope to be able to apply these 
results to other classes of systems or more complex ones, e.g. with delays. Complexity reduction is also an important 
and interesting research direction. 
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